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Abstract 

An  experiment  with  a  simulated  microeconomic  system  demonstrates  that  the  decision- 
making processes  of  human  subjects  can  produce  deterministic  chaos.   Participants  managed 
a  commodity  production-distribution  system  to  minimize  costs.   Performance,  however,  was 
systematically  suboptimal.   A  model  of  the  subjects'  decision  rule  is  proposed.   Econometric 
estimates  show  the  model  is  an  excellent  representation  of  the  actual  decisions.    Simulation 
of  the  estimated  rules  yields  stable,  periodic,  quasiperiodic,  and  chaotic  solutions.   Analysis 
of  the  parameter  space  reveals  a  complex  structure  including  mode-locking,  a  devil's  stair- 
case, and  fractal  basin  boundaries.   Implications  for  modeling  human  systems  are  explored. 
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The  prevalence  of  deterministic  chaos  and  other  nonlinear  phenomena  in  physical, 
chemical,  and  biological  systems  has  prompted  speculation  that  these  dynamics  may  occur  in 
economic,  managerial,  and  other  human  systems  as  well  (1).   Indeed,  numerous  models  have 
shown  how  economic  systems  might  produce  chaos  (2).   But  the  significance  of  these 
important  theoretical  developments  hinges  on  whether  the  chaotic  regimes  lie  in  the  realistic 
region  of  parameter  space.   Further,  the  decision  rules  postulated  in  these  models  have  not 
been  tested  experimentally.  Despite  intriguing  efforts  to  identify  chaos  in  economic  time 
series  (3),  it  is  difficult  to  resolve  such  issues  by  empirical  means  alone  (4).  Data  series  are 
often  unavailable,  or  too  short  relative  to  the  frequencies  of  interest.   Measurement  error  and 
process  noise  in  social  and  economic  data  are  significant,  further  complicating  empirical 
analysis  (5).   A  complementary  approach  is  based  on  laboratory  experiments  in  which  models 
provide  a  simulated  environment  for  the  study  of  decision-making  (6).  We  report  here  the 
results  of  one  such  experiment  which  demonstrates  that  the  decision-making  processes  of 
human  subjects  can  produce  deterministic  chaos  in  a  common  economic  context. 

The  experiment  creates  a  simulated  industrial  production-distribution  system.    Such 
systems  give  firms  flexibility  by  establishing  a  decentralized  network  of  inventories  to  buffer 
unanticipated  variations  between  the  demand  for  and  production  of  commodities  and 
manufactured  goods.    Production-distribution  systems  have  long  been  associated  with 
business  cycles  and  other  fluctuations  in  industrial  economies  (7).   The  experiment  here,  the 
"Beer  Distribution  Game,"  is  a  role-playing  simulation  of  a  typical  production-distribution 
system.   Developed  at  MIT  to  introduce  students  of  management  to  economic  dynamics  and 
computer  simulation,  the  game  has  been  widely  used  for  nearly  three  decades  (8). 

The  experiment  is  conducted  on  a  board  which  portrays  in  a  simplified  fashion  the 
production  and  distribution  of  beer  (figure  1).   Beer  is  represented  by  markers  which  are 
manipulated  by  the  players  as  the  game  proceeds.   Each  brewery  consists  of  four  sectors: 
retailer,  wholesaler,  distributor,  and  factory  (R,  W,  D,  F).   One  subject  manages  each  sector. 
Customer  demand  is  represented  by  a  deck  of  cards.   Each  simulated  week  customers  order 
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beer  from  the  retailer,  who  ships  the  beer  requested  out  of  inventory.  The  retailer  in  turn 
orders  from  the  wholesaler,  who  likewise  ships  out  of  inventory.   Similarly,  the  wholesaler 
orders  and  receives  beer  from  the  distributor,  who  in  turn  orders  and  receives  beer  from  the 
factory.  The  factory  brews  the  beer.  At  each  stage  there  are  shipping  and  order  receiving 
delays.   Each  week  the  subjects  must  decide  how  much  to  order  from  their  immediate  supplier 
(factories  decide  how  much  beer  to  brew). 

Subjects  seek  to  minimize  cumulative  costs  incurred  during  the  game  (36  simulated 
weeks).   Inventory  holding  costs  for  each  sector  are  $.50/case/week,  and  backlogs  of  unfilled 
orders  cost  $  1 .00/case/week.   Subjects  must  therefore  keep  their  inventory  low  while 
avoiding  backlogs.  Due  to  the  order  receiving  and  shipping  lags,  a  substantial  supply  line  of 
unfilled  orders  may  build  up.  The  lag  in  receiving  beer  may  vary:  if  the  wholesaler  can  cover 
the  retailer's  orders,  the  retailer  will  receive  the  beer  desired  after  four  weeks.   But  if  the 
wholesaler  has  run  out,  the  retailer  must  wait  until  the  wholesaler  receives  additional  beer. 
Only  the  factory,  the  primary  producer,  faces  a  constant  delay  in  acquiring  inventory  (for 
simplicity  there  is  no  limit  to  production  capacity).   Although  the  experimental  system  is 
highly  simplified,  it  nevertheless  shares  many  features  with  real  economies,  including 
multiple  feedbacks,  nonlinearities,  and  time  lags.   The  nonlinearities  arise  from  non- 
negativity  constraints  on  orders  and  shipments:   shipments  normally  equal  incoming  orders, 
but  if  inventory  is  depleted,  shipments  must  equal  zero,  and  the  unfilled  orders  are  held  in  the 
backlog  for  future  delivery. 

The  experiment  follows  standard  protocols  (6,  9)  and  is  detailed  in  (10).  Each  trial  is 
initialized  in  equilibrium.  Each  inventory  contains  12  cases.  Customer  demand  is  initially  4 
cases  per  week.   The  system  is  disturbed  by  an  unannounced  step  increase  in  customer 
demand  to  8  cases  per  week  in  week  5. 

Forty-eight  trials  (192  subjects)  were  analyzed.  Trials  in  which  subjects  made 
accounting  errors  were  discarded,  leaving  eleven  (44  subjects).  This  subset  generated 
slightly  lower  costs  than  the  full  sample,  indicating  that  these  subjects  understood  and 
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performed  best  in  the  experiment.    Participants  were  graduate  and  undergraduate  students  at 
MIT  and  senior  executives  from  a  number  of  U.S.  firms. 

Results  are  summarized  in  table  1;  figure  2  shows  a  typical  trial.   The  participants' 
behavior  is  significantly  suboptimal.  Given  the  complexity  of  the  system  it  is  not  surprising 
that  subjects  are  unable  to  minimize  costs.    More  interesting,  the  results  exhibit  strong 
regularities,  suggesting  subjects  used  a  common  heuristic  in  ordering.   1.   Oscillation:  Orders 
and  inventory  are  dominated  by  a  large  amplitude  fluctuation.  On  average  21  weeks  are 
required  to  recover  initial  inventory  levels.  2.  Amplification:  The  amplitude  and  variance  of 
orders  rises  steadily  from  customer  to  retailer  to  factory.   Customer  orders  increase  from  4  to 
8  cases/week;  responding  to  this  disturbance,  factory  orders  rise  from  4  to  a  peak  averaging 
32  cases.  3.  Phase  lag:  Orders  tend  to  peak  later  as  one  moves  from  retailer  to  factory  (not 
surprising  since  orders  propagate  upstream  through  decision-making  and  processing  delays). 

We  next  estimate  a  model  of  the  subjects'  decision  rule  (see  (10)  for  details).   We 
hypothesize  that  subjects  utilize  the  anchoring  and  adjustment  heuristic  to  determine  orders 
(11).  Anchoring  and  adjustment  is  a  common  judgmental  strategy  in  which  an  unknown 
quantity  is  estimated  by  first  recalling  a  known  reference  point  (the  anchor)  and  then 
adjusting  for  factors  which  may  be  less  salient  or  whose  effects  are  obscure.   For  example,  a 
person's  weight  may  be  estimated  by  recalling  the  subject's  own  weight  and  adjusting  for 
differences  in  height,  girth,  etc.  Anchoring  and  adjustment  has  been  documented  in  a  wide 
variety  of  decision-making  tasks  (11,  12). 

Here  the  anchor  is  expected  demand  D^,  the  subject's  forecast  of  incoming  orders. 
Replacement  of  incoming  orders,  ceteris  paribus,  maintains  inventory  at  its  current  value. 
The  asymmetric  cost  function  means  subjects  should  maintain  their  inventory  at  a 
nonnegative  value  S'.   Thus  an  adjustment  to  the  anchor  arises  from  the  discrepancy  between 
the  desired  stock  S'  and  actual  stock  S.   Subjects  may  also  consider  a  fraction  P  of  the  supply 
line  of  unfilled  orders.  The  supply  line  SL  is  the  accumulation  of  orders  placed  but  not  yet 
received.   Finally,  orders  O  must  be  nonnegative  (order  cancellations  are  not  allowed).   Thus, 
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Ot=  MAX(0,  D^t  +  a(S'  -  Sf  pSLj))  (1) 

where  the  stock  adjustment  parameter  a  represents  the  fraction  of  any  discrepancy  corrected 
each  time  period. 

The  adjustment  term  a(S'  -  Sf  pSLj)  creates  two  negative  feedback  loops  which    " 

regulate  the  stock.   Discrepancies  between  the  desired  and  actual  stock  induce  additional 
orders  until  inventory  reaches  the  desired  value.  Likewise  an  insufficient  quantity  of  goods 
on  order  induces  additional  orders  until  the  pipeline  is  filled  -  depending  on  p.  If  p  =  0,  then 
orders  placed  are  forgotten  until  they  arrive,  causing  overordering  and  instability.    If  P  =  1, 
then  subjects  fully  recognize  the  supply  line  and  do  not  double  order.  While  p=l  is  optimal,  a 
prior  experiment  showed  P«l  for  many  subjects  in  a  similar  task  (13). 

Expected  demand  EK  represents  each  participant's  forecast  of  incoming  orders. 
Adaptive  expectations  are  assumed: 

D^  =  eDt.i  +  (i-e)D^t-i,   o<e<i.  (2) 

Consistent  with  behavioral  decision  theory  (11,  12,  14),  the  rule  utilizes  information 
locally  available  to  the  decision  maker  and  does  not  presume  managers  have  global  knowl- 
edge of  the  system  structure  or  the  cognitive  capability  to  solve  for  optimum  performance. 

The  parameters  (0,  a,  P,  and  S')  were  estimated  for  each  participant  by  nonlinear 
least  squares  subject  to  the  constraints  0<6<1  and  a,  P,  S'  >0  (table  2).   The  explanatory 
power  of  the  proposed  rule  is  excellent:  the  mean  R^  is  71%;  R2  is  less  than  50%  for  only  6  of 
44  subjects.   A  large  majority  of  the  estimated  parameters  are  significant,  and  the  estimated 
parameters  are  systematically  related  to  performance  (10). 

Though  the  subjects'  disequilibrium  response  is  suboptimal,  we  expected  that  their 
decision  rules  would  be  stable  and  swiftly  return  the  system  to  a  low-cost  equilibrium.   To 
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test  this  hypothesis  we  simulated  the  experimental  system  using  each  set  of  estimated 
parameters  (15).   Thirty  parameter  sets  (68%)  do  indeed  produce  stable  behavior.   However, 
one  periodic  and  three  quasiperiodic  solutions  appear.  Ten  (23%)  yield  aperiodic  (chaotic) 
behavior.   Figure  3  shows  several  characteristic  simulations.   Compare  the  shape  of  the 
strange  attractor  produced  by  the  parameters  of  subject  4  (0,  a,  P,  S'  =  1.0,  0.65,  0.40,  15)  to 
subject  36  (0,  a,  (3,  S'  =  0,  0.3,  0.15,  17).  Note  the  long  and  apparently  chaotic  transient 
produced  by  the  (stable)  parameters  for  subject  27  (0,  a,  p,  S'  =  0.2,  0.3,  0.05,  8). 

To  explore  the  structure  of  the  parameter  space,  we  set  0  and  S'  at  representative 
values  of  .25  and  17,  respectively,  and  varied  a  and  P  over  the  managerially  meaningful 
interval  [0,1].   The  space  (figure  4)  is  complex,  including  stable,  periodic,  quasiperiodic,  and 
chaotic  solutions.    The  unstable  solutions  arise  from  the  nonlinear  coupling  of  several  oscil- 
latory feedbacks  created  by  the  multiple  inventories  and  time  delays  in  the  system.   These 
oscillators  are  coupled  through  the  availability  of  inventory  in  the  distribution  chain.  High 
frequencies  are  produced  when  there  is  sufficient  inventory  so  that  orders  can  be  filled  by 
each  sector's  immediate  supplier,  e.g.,  when  the  retailer's  orders  aie  filled  out  of  the  whole- 
saler's inventory.   Low  frequencies  arise  when  inventories  are  inadequate,  forcing  customers 
to  wait  for  the  factory  to  receive  orders,  produce  the  beer,  and  finally  ship  it  downstream. 
Periodic  solutions  arise  when  the  parameters  of  the  heuristic  are  such  that  the  ratio  of  these 
frequencies  is  rational;  quasiperiodic  solutions  arise  when  this  winding  number  is  irrational, 
and  chaotic  solutions  arise  when  the  individual  oscillators  are  sufficiently  unstable  (16). 

In  general,  larger  values  of  a  and  smaller  values  of  p  are  destabilizing  (figure  4).  To 
the  extent  a  subject  ignores  the  supply  line  of  unfilled  orders  (P<1),  a  stock  shortfall  causes 
orders  to  be  placed  each  period  even  after  sufficient  orders  are  in  the  supply  line,  leading  to 
excess  inventories  and  oscillation.   Such  overordering  is  exacerbated  by  more  aggressive 
stock  adjustment  (a>0.4)  since  more  will  be  ordered  in  response  to  a  given  stock  shortfall. 
However,  the  transitions  between  modes  are  not  monotonic.   Note  particularly  the  two 
narrow  fjords  of  stable  solutions  which  snake  between  the  periodic  and  chaotic  solutions. 
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Figure  5  magnifies  a  region  of  parameter  space  at  the  transition  from  stable  to 
unstable  solutions.   The  alternating  bands  of  periodic  and  quasiperiodic  solutions  arise  from 
mode-locking  between  the  various  oscillators  and  form  a  devil's  staircase  (16).   Further 
magnification  reveals  the  characteristic  fractal  steps  of  the  staircase.   The  structure  of  the 
staircase  is  complicated;  however,  by  the  fingers  of  stable  solutions  which  cut  across  the 
staircase.   These  regions  of  stability  appear  to  have  fractal  boundaries  as  well.   A  transition 
directly  from  stability  to  chaos  also  exists  and  also  appears  to  be  fractal.   Figure  6  varies  P 
for  a=0.90.  As  P  decreases  from  .444  to  .434  the  system  moves  irregularly  from  stability  to 
chaos.   Successive  magnification  reveals  the  self-similarity  of  the  fractal  boundary. 

It  is  common  in  the  social  sciences  to  assume  that  decision-making  behavior  and  thus 
the  dynamics  of  human  systems  are,  if  not  optimal,  then  at  least  stable.   These  results  show 
that  formal  rules  which  characterize  actual  managerial  decision  making  can  produce  an 
extraordinary  range  of  disequilibrium  dynamics,  including  chaos.   Such  complexity  raises  a 
number  of  issues  for  social  scientists.  Policy  interventions  often  imply  changes  in  the 
parameters  of  a  decision  rule  or  model.  How  can  policy  analysis  be  conducted  if  the  "policy 
space"  contains  fractal  basin  boundaries?   In  such  systems  changes  on  the  margin  may 
produce  unpredictable  qualitative  changes  in  behavior.   Experience  may  not  transfer  to 
circumstances  which  differ  only  slightly.  Do  robust  principles  of  policy  design  exist  in  such 
systems? 

Even  though  deterministic  cause-effect  relations  exist  in  chaotic  systems,  it  is 
impossible  to  predict  the  effects  of  small  changes  in  initial  conditions  or  parameters.  Does 
such  'randomness'  slow  the  discovery  of  cause  and  effect  by  agents  in  the  economy  and  thus 
hinder  learning  or  evolution  towards  efficiency?  Indeed,  does  learning  alter  the  parameters  of 
decision  rules  so  that  systems  evolve  towards  or  away  from  the  chaotic  regime?  Much 
further  development  of  theory  and  experiment  are  required  to  answer  these  questions  and 
assess  the  significance  of  chaos  in  human  systems. 
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Figure  1.   "Beer  Distribution  Game"  board.    Subjects  carry  out  five  steps  each  simulated 
week: 

1.  Receive  inventory  and  advance  shipping  delays. 

2.  Fill  orders. 

3.  Record  inventory  or  backlog. 

4.  Advance  the  order  slips. 

5.  Place  orders. 

Subjects  make  their  decisions  in  step  5;  steps  1-4  handle  bookkeeping  and  other  mechanics. 


Figure  2.  Typical  experimental  results.  Top:  customer  orders  increase  from  4  to  8 
cases/week  in  week  5.   Middle:  The  resulting  orders  placed  by  subjects  (from  bottom  to  top, 
Retailer,  Wholesaler,  Distributor,  Factory).   Bottom:  Effective  inventory  levels  (Effective 
inventory=Inventory-Backlog).   Tick  marks  on  y-axes  denote  10  units.   Note  the  oscillation, 
amplification,  and  phase  lag  as  the  disturbance  propagates  from  customer  to  factory. 


Figure  3.   Simulation  of  the  decision  rule  with  estimated  parameters.   Top:  phase  portraits 
showing  distributor  inventory  vs.  factory  inventory.   Weeks  1000-1600  shown  to  remove 
transient;  flow  is  clockwise.  Bottom:  The  parameters  for  subject  27  produce  a  long  chaotic 
transient  (Retailer  orders  shown). 


Figure  4.     Parameter  space  of  the  simulated  system.   The  stock  and  supply  line  adjustment 
parameters  a  and  P  are  varied  over  the  interval  [0,1]  in  increments  of  .025,  holding  0=.25  and 
S'=17.   Green=stable;  red=periodic;  blue=quasiperiodic;  orange=chaotic.  Note  the  two  fjords 
of  stable  solutions  separating  regions  of  unstable  behavior. 


Figure  5.  Magnification  of  parameter  space  by  a  factor  of  100  compared  to  figure  4. 
Green=stable;  red=periodic;  blue=quasiperiodic;  yellow=period-doubling.  The  periodic 
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solutions  have  rational  winding  number,  indicating  the  high  frequency  oscillator  locks  into  a 
rational  multiple  of  the  low  frequency  cycle.  The  numbers  3,  4,  5,  etc.  indicate  the  number  of 
low  frequency  cycles  completed  before  the  system  repeats. 


Figure  6.  Transition  from  stability  to  chaos  as  P  varies,  holding  a=.90.   Green=stable; 
orange=chaotic.   Successive  magnification  reveals  the  irregularity  characteristic  of  a  fractal 
boundary.  Rightmost  column  of  figure  reflects  magnification  of  10^  compared  to  figure  4. 
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Table  1.   Summary  of  experimental  results.   Data  are  means  of  11  trials.   Actual  costs  exceed 
optimal  levels  by  a  factor  of  10,  a  highly  significant  margin. 

Total  Retailer       Wholesaler      Distributor     Factory 


COSTS 
Mean 

Benchmark* 
Ratio 

t-statistic: 

Hq:  Mean  cost  = 

Benchmark 


$2028 

$383 

$635 

$204 

$46 

$50 

9.9 

8.3 

12.7 

$630 

$380 

$54 

$54 

11.7 

7 

8.7  4.9  5.9  6.9  9.7 

p<.000+       p<.001  p<.000+  p<.000+      p<.000+ 


Customer      Retailer       Wholesaler      Distributor     Factory 


TIMING 
Time  to  recover  initial  inventory  (weeks)      N/A 


AMPLIFICATION 

Peak  Order  Rate  (cases/week) 
Variance  of  Order  Rate  (cases/week)^ 

PHASE  LAG 
Date  of  Peak  Order  Rate  (week) 


24 


16 


23 


16 


22 


21 


16 


8 

15 

19 

27 

32 

1.6 

13 

23 

45 

72 

20 


'Benchmark  costs  are  the  minimum  costs  produced  by  simulation  of  the  proposed  decision  rule  for 
orders  (eq.  1-3)  and  are  an  upper  bound  estimate  of  optimal  performance  in  the  experiment. 
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Table  2.  Summary  of  estimation  results  for  44  subjects.  The  ordering  heuristic  is  given  by  eqs. 
1-2.   0  determines  the  speed  of  adjustment  of  the  demand  forecast;  a  determines  the 
response  to  inventory  shortfalls;  (3  is  the  fraction  of  the  supply  line  accounted  for  by  the 
subjects;  S'  is  the  desired  stock.   See  (10)  for  details. 

Parameter:  6  a  P  S'  R^ 

Mean  (std.  dev.)  0.36(0.35)  0.26(0.18)  0.34(0.31)  17  (9)         0.71    (0.22) 

Median:  0.25  0.28  0.30  15                 0.76 

Minimum:  0.00  0.00  0.00  0                 0.10 

Maximum:  1.00  0.80  1.05  38                0.98 
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Appendix:  Equations  for  Production-Distribution  System 

The  equations  are  presented  in  the  sequence  of  calculation  used  in  the  experiment  and 
simulations  (figure  1).  Initial  equilibrium  conditions  are:  all  inventories  =  12  units,  all 
backlogs  =  0,  all  other  variables  =  4  units  (orders  placed,  incoming  orders,  production  delays, 
shipping  delays,  and  expected  orders).   Customer  demand  is  initially  4  and  rises  to  8  in  week 
5.    Subscripts  denote  the  sector  of  the  distribution  chain  (Retailer  R,  Wholesaler  W, 
Distributor  D,  and  Factory  F).   The  computer  programs  for  the  simulations  and  parameter 
estimation  are  available  from  the  first  author. 

Step  1:  The  contents  of  the  shipping  delay  (Dl)  immediately  to  the  right  of  each 
sector's  inventory  (I)  are  added  to  inventory.   The  contents  of  the  shipping  delay  on  the  far 
right  (Dl)  are  moved  into  Dl.  Factories  advance  the  production  delays  Dip  and  D2p  in  the 
same  fashion. 

Ij(t)  =  Ij(t-l)  +  Dlj(t-l)  (Al) 

Dlj(t)  =  D2j(t-l)  (A2) 

Step  2:   Retailers  examine  the  top  card  on  the  Customer  Order  deck  (CO);  all  others 
examine  Incoming  Orders  (10).   All  sectors  fill  orders.   The  Shipment  rate  S  must  equal  the 
new  orders  plus  any  Backlog  B  from  the  prior  period,  to  the  extent  inventory  permits.  The 
retailer's  shipments  go  directly  to  the  customer  and  leave  the  system.   Shipments  of  all 
others  are  placed  in  the  shipping  delay  D2  of  the  downstream  sector  (A4).    Shipments  also 
reduce  inventory  (A5). 

SR(t)  =  MIN(CO(t)  -f  BR(t-l),  lR(t))  (A3) 

Si(t)  =  MIN(IOi.i(t-l)  +  B|(t-l),Ii(t)),      1  =  W,D,F  (A3') 

D2k(t)  =  Sk^i(t),       k  =  R,W,D  (A4) 

Ij(t)  =  Ij(t)-Sj(t)  •  (A5) 

Step  3:   All  sectors  record  inventory  or  backlog.  The  net  change  in  backlog  is  the  differ- 
ence between  incoming  orders  and  shipments.   Effective  Inventory  EI  is  inventory  less  backlog. 
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BrCD  =  BR(t- 1 )  +  CO(t)  -  SrCD  (A6) 

Bi(t)  =  B,(t-l)  +  IO,.i(t-l)-  S,(t),      1  =  W,D,F  (A6') 

EIj(t)  =  Ij(t)-  Bj(t)  (A7) 

Step  4:   Each  sector  advances  the  slip  containing  last  week's  Order  O  to  incoming 
orders.  The  factory  puts  last  week's  order  in  the  top  production  delay  D2p. 

I0,,(t)  =  0j,(t-1),       k  =  R,W,D  (A8) 

D2F(t)  =  OF(t-l)  (A8') 

Step  5:  Each  sector  places  orders.   First  the  Supply  Lines  SL  are  calculated.   The 
supply  line  is  the  sum  of  units  in  the  two  shipping  delays,  the  backlog  of  the  supplier  (if  any) 
and  orders  placed  the  previous  week.   Since  the  factory  is  the  primary  producer,  its  supply 
line  is  simply  the  contents  of  the  production  delays. 

SLk(t)  =  Dli.(t)  +  D2i,(t)  +  Bk^i(t)  +  I0j,(t),       k  =  R,W,D  (A9) 

SLF(t)  =  DlF{t)  +  D2F(t)  (A9') 

Each  sector's  forecast  of  incoming  orders  (expected  Demand)  D^  is  formed  by  adaptive 
expectations,  with  smoothing  parameter  6. 

D-(t)  =  ejlOj(t- 1 )  +  ( l-ej)D.(t- 1 )  (AlO) 

Finally,  orders  are  the  given  by  demand  forecast  adjusted  by  a  fraction  a  of  the  difference 
between  the  desired  Stock  S'  and  the  effective  inventory,  including  a  fraction  p  of  the  supply 
line.     Orders  must  be  nonnegative. 

Oj(t)  =  MAX(0,D-(t)  +  aj(S'j-  Elj(t)-  PjSL/t)))  (All) 
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Trial  &  Position 


Appendix.  Estimated  parameters,  Beer  Distribution  Game 

e  a  P  S'  r2 


Mode 


10 


11 


R 
W 
D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 

R 

W 

D 

F 


0.90 

0.00 

0.15 

1.00 

a 

0.00 

0.00 

0.00 

0.25 

c 

0.00 

0.00 

0.90 

a 

0.25 

0.40 

0.30 

0.05 

0.25 

0.10 

0.95 

a 

0.20 

b 

0.25 

0.05 

0.30 

0.15 

0.55 

a 

0.95 

0.50 

a 

0.20 

a 

0.80 

b 

0.50 

0.40 

a 

1.00 

a 

0.55 

a 

0.05 

0.20 

0.30 

a 

0.00 

1.00 

0.05 

0.15 

0.40 

a 

0.75 

0.00 

0.05 

0.95 

a 

0.10 

0.25 

a 

0.05 

a 

0.65 

a 

0.40 

a 

0.40 

a 

0.30 

a 

0.25 

a 

0.20 

a 

0.15 

a 

0.30 

a 

0.30 

a 

0.35 

a 

0.05 

a 

0.35 

a 

0.25 

a 

0.35 

a 

0.15 

a 

0.20 

a 

0.35 

a 

0.30 

a 

0.20 

a 

0,05 

0.65 

a 

0.15 

a 

0.00 

0.30 

a 

0.00 

0.05 

0.10 

a 

0.15 

a 

0.80 

a 

0.30 

a 

0.00 

0.10 

a 

0.30 

a 

0.00 

0.30 

a 

0.60 

a 

0.35 

a 

0.35 

a 

0.25 

a 

0.50 

a 

0.30 

b 

0.65 

a 

20 

a 

0.20 

0.50 

a 

27 

a 

0.86 

0.35 

14 

0.74 

0.40 

a 

15 

a 

0.84 

0.10 

a 

7 

a 

0.67 

0.75 

a 

30 

a 

0.92 

0.10 

a 

10 

a 

0.88 

0.10 

9 

a 

0.87 

0.00 

25 

a 

0.57 

0.50 

a 

38 

a 

0.11 

0.20 

a 

10 

a 

0.61 

0.00 

18 

a 

0.73 

0.45 

a 

15 

a 

0.43 

0.00 

30 

c 

0.76 

1.00 

a 

18 

a 

0.86 

0.00 

19 

a 

0.89 

0.65 

a 

13 

a 

0.60 

0.55 

a 

14 

a 

0.79 

0.30 

a 

19 

a 

0.94 

0.55 

a 

24 

a 

0.73 

0.65 

a 

31 

a 

0.58 

0.35 

a 

27 

a 

0.82 

0.25 

15 

0.32 

0.00 

9 

a 

0.75 

0.00 

9 

a 

0.75 

N/D 

N/D 

0,87 

0.05 

a 

8 

a 

0.98 

N/D 

N/D 

0.60 

0.60 

6 

0.10 

0.30 

a 

16 

a 

0.81 

0.80 

a 

14 

a 

0.73 

0.00 

9 

a 

0.87 

0.45 

a 

5 

a 

0.89 

N/D 

N/D 

0.23 

0.90 

a 

12 

a 

0.94 

0.15 

c 

17 

a 

0.87 

WD 

N/D 

0.76 

0.20 

a 

20 

a 

0.76 

0.35 

a 

0 

0.69 

1.05 

a 

32 

a 

0.95 

0.00 

4 

a 

0.83 

0.05 

18 

a 

0.72 

0.00 

15 

a 

0.84 

0.20 

26 

a 

0.66 

s 

PI 

s 

c 

c 

s 

c 

s 

s 

s 

Q 

s 
s 
s 
s 
s 

Q 
S 
S 
S 

s 
s 
s 
c 
s 
s 
s 
s 
s 
s 
s 
c 

Q 

s 
s 
c 
s 
s 
c 
s 
c 
c 
c 
s 


N/D:  Not  Defined 

Significant  at  a:  .005;  b:  .01 ;  c:  .025  level  (1 -tailed  t-test  [since  parameters  must  be  >  0]) 

Mode:  S=Stable,  Pn=Limit  cycle  of  period  n,  Q=Quasiperiodic,  C=Chaotic 
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